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Objectives 

The objective of this paper is to study the stability characteristics of a 
bearing tester, tfe verify our conclusions using nuaerical simulations of a 
realistic model. 
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Abstract. He study the properties of the solutions of a system of four 


coupled nonlinear differential equations that model the behavior of the 
rotating shaft of a bearing tester. In particular, we show how bounds for the 
solutions of these equation can be obtained from bounds for the solutions of 
the linearized equations. By studying the behavior of the Fourier transforms 
of the solutions, we are also able to predict the approach to the stability 
boundary. These conclusions are verified by means of numerical solutions of 
the equations, and of power spectrum density (PSD) plots. 
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1 . Introduction, 



In this study we continue the investigation of the properties of the 
solutions of mathematical models of rotating aachinery initiated by Day [1]. 
Both Day and this author [2] have studied the behavior of a simple Jeffcott 
aodel with deadband, viz. a system of coupled differential equations that 
represent the behavior of a rotating shaft. 

The purpose of this paper is to examine the properties of the solutions 
of a model of a bearing tester. This is a device designed to estimate the 
life expectancy of bearings under realistic conditions of loads and 
acceleration in cryogenic fluids. Our study will help determine safety 
margins for its operation. 

We consider a bearing tester with two seals and two bearings with 
deadband. A sketch of this mechanical system can be seen in Fig. 1. 

2. General Theory 

2.1 Derivation of the Bearing Tester Equations. 

We assume that the shaft is rotating with angular velocity u along an 
axis close to the x-axls, that both bearings are at the same distance a from 
the center of symmetry of the shaft, that both seals are at the same distance 
b from this center of symmetry, and that the shaft cannot move in the 
direction of the x-axis. We also assume that both seals have the same damping 


C , stiffness K , and cross coupling stiffness Q , and that both bearings 
8 8 8 

have the same stiffness K^. 

For j =■ 1, 2, let 4^ denote the magnitude of the deadband at bearing 
1; let v . and v . describe the displacement of the center of the shaft at 

y i z j 

bearing 1, and let w . and w . be similarly defined for the seals (see 

y 3 z j 

2 2 1/2 

Fig. 2). Let a denote the mass of the shaft. If r ■ (v + v ) , 

3 73 Z J 
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h^(t) » 1 If r^ < 6^, and h^(t) « 6^/r^ if r^ > 6^ , then the equations 
that describe the aoveaent of the shaft are the following : 


V H l (tll ',l + K b“' h 2 <t)1 V + + V 1 + V“zl * “z2 ] + 


c s i»;i + y + (,/2,lv ;i + ' g i (t> 


a) 


K b 11 - h l <t)1 ’z 1 + K b ll - h 2 (t »''z2 + V“zl + “z2« * V"yl * V 1 + 

c JKi * “k 1 + <* /2)lv ^i + v "z 2 ) - s 2 (t) 


( 2 ) 


+bC s [^2 * «’! 


+ aK b U-h 2 (t)) V + K s [» y2 - „ yl ] + bQjw^ 
] + <I 2 /2a)[.- 2 - y- yl ] - (»I 1 /2a)[,; 2 - y^l 


“zl 1 


■j(t) 


(3) 


-aK b [l-yt)l, zl * aK b U-b 2 (t)]v z2 + <*„l « z2 - » zl l ' *VV * V 1 

(4) 

+ “sKl - “il 1 + «2 /2 " >l ’^2 - ’zl 1 + ( “ I l /2 *) ( ’i2 * ’il 1 * *2<‘> 


V , + ▼ n * W 2 + W o» 

yi y2 yi y2 


v. - v „ =» w _ - w, 
yl y2 yl y2 


V + V „ = w , + w „ 
zl z2 zl z2 


v , - v - w , - w 
zl z2 zl z2 


(5) 

( 6 ) 


where I is the axial inertia, I is the rotational inertia about the axis 
1 2 

transversal to the shaft, and g^Ct), g^t), a^(t), a^t) are given as 
follows: 

g]L (t) - u 2 «I(e yl + e y2 )cos tot - (e zl + e z2 > sin totl 

g^t) - 0) «l( e zl + e z 2^ cos ut + ^ e yl + e y2^ Si “ 

■jCt) - u) 2 al(e yl - e y2 ) cos tot - (e zl - e z2 ) sin tot) 

^(t) - «i) 2 «l(e zl - e z2 ) cos tot + (e yl - e y2 > sin tot]. 
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where e - e _ , e , and e _ represent the mass Imbalance. In our analyses 
yl y2 zl z2 

we shall assume that g^, m^, and m^ are arbitrary continuous and 

bounded functions. 

Note that (1) and (2) are force equations, (3) and (4) are moment 
equations, and (5) and (6) are derived from the symmetry assumptions on 


bearings and seals. Setting 
g(t) - g 1 (t) + ig 2 (t) and m(t) 


v . + iv ., w. « w + iw , 

y i z j j y j z j 

m^(t) + im 2 (t), we obtain: 


K b [1 - h l ( t)'» 1 + K b H-» 2 (t)W 2 ♦ (K„ - «)„>[»! + « 2 1 
+ C g [w' + + (■/2)Ivj + v;] =* g(t). 

-a^d-hjCtnVj + aK b U-l, 2 (t)]v 2 + b(K a - lQ 8 )[- 2 - + bC.I.J 

+ (I 2 /2a)[v 2 - vjl ♦ - vjl = m(t). 


( 7 ) 


w'l 

1 J 

( 8 ) 


and 


V 1 + V 2 


W 1 + 


2 * 



(9) 


In view of (9), if we set v « v + v^, u * v 2 - v^, q^ * h^v^, and 
q 2 » t^v^ (7) and (8) can be written in the following form: 


(m/2)v" + C v* + (K. + K - iQ )v - R. [q, + q_] - g(t), 
s b 8 8 b I 2 


and 


(l /2a)u" + (bC + itiil /2a)u* + (aK + bK - ibQ )u + aK [q - q ] * m(t). 
2 Sl D88 DlZ 

Thus, setting * 2C g /m, C 2 * (2abC g + iti>I^)/I 2 ,K^ * 2K^/m, K 2 * (2a K^)/^, 
\ x =* 2K g /m, A 2 - 2abK g /I 2 , B x - 2Q g /m, B 2 - 2abQ g /I 2 , M x “ A x + \ “ iB x » 
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( 10 ) 


M 2 - A 2 + K 2 - IB f t (t) « (2/a) g(t), f 2 (t) = (2a/I 2 )a(t), and 

P 1 * ’l + V P 2 * "l ' ’2’ 

we finally obtain the Bearing Tester equations: 

v" + (^v* + MjV - K x p x - f^t), (11) 

and 


u" + C 2 u* + M 2 u + K 2 p 2 - f £ (t). (12) 

We shall as suae that B , B 2 , C^, C 2 , K^, & 2 are positive, and A^ t 

A 2 and t are nonnegative. 

Since Vj - (l/2)(v - u) and v 2 - (l/2)(v + u), p^t) and p 2 (t) can 
be expressed in teras of v and u using (10) and the following 
representations for q^(t) and q 2 (t): 

(v-u)/2, if I v— u| < 26, 

q.(t) ■ » (13) 

« 1 (v ta ni)/ jv-u| , if |v-u| > 

and 


(v+U>/2, if |v*u| < 2« 2 

6 2 (v+u)/ jv+uj , If J v+u j > 26 2 « 


(14) 


2.2 Existence, uniqueness, and representation foraulas. 


We have transformed the system of equations (1) - (6) into the equivalent 
system (11), (12). This is a system of coupled nonlinear differential 
equations similar to the Jeffcott equations we studied in [2]. The existence 
and uniqueness of their solutions (and therefore, of the solutions of the 
original system), follow by the same argument employed for the Jeffcott 
equations, and need not be repeated here. 


Let Qj, * -4( A j + Kj), 

*1 ■ »‘ w «, * «£ - K> 1/2 i ,/2 . 


(15) 
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( 16 ) 

(17) 


“l * “ C l 1/2 » 

\ ma l* IBj, 


+ CjJ/2, 

“i ’ 1B 1- 


Then, as In [2] , 1c Is readily seen that A^ and X^ are the solutions of 


the characteristic equation X + C^A + = 0, and therefore 


v h - Cj exp(X 1 t) + c 2 exp(X 2 t) (18) 

Is the general solution of 

v" + CjV' + M x v - 0, (19) 

Similarly, If Yj * <* 2 + and ^2 " a 2 + 1B 2 are the 8 °l utl<m8 the 

2 

characteristic equation y + C^r + M,, = 0, it Is clear that 

u fc « dj exp(Yjt) + d 2 exp(r 2 t) (20) 

Is the general solution of 

u" + C 2 u' + l^u - 0. (21) 

Without loss of generality, we shall always assume that a^. If C 2 Is 

real (l.e. , If oil^ * 0), then formulas similar to (IS), (16) and (17) obtain 
for and y 2 - 

If Vp and Up are particular solutions of the linearllzed Bearing 
Tester equations 

v" + Cj v' + MjV - f (t) (22) 

and 

u" + C 2 u' + M 2 u - f 2 (t), (23) 

then, setting v t * v h + v p » u t “ u h + U p* 

G 1 (t) - (Xj - X 2 ) 1 [exp(X 1 t) - exp(X 2 t)], 

G 2 (t) - (Yj “ Y 2 ) 1 lexp(Y 1 t) - exp(Y 2 t)l, 
and proceeding as in [2], we readily deduce that (11) and (12) are equivalent 
to the following nonlinear Vol terra Integral equations of convolution type: 
v(t) - v t (t) + P x (t) (24) 

and 
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u(t) * U^(t) + P 2 (t), 


(25) 


where the perturbation tens ( t) are given by: 


P j (t> " K j f 0 Gjtfc-’OPj <*>**» j “ 1. 2. (26) 

Thus, the existence and uniqueness of the solutions of (11) and (12) also 
follow froa the existence and uniqueness of the solutions of (24) and (25) 
(cf. e.g. [3], [4]). 


2.2 Bounds 


Let D^(t) * f^jGj(t-x) jdx, and let 6 = + $ 2 . Since (10), (11), (13) 

and (14) imply that 

|Pj(t)| < «, j - 1, 2, (27) 

we readily conclude that 

|P j (t)| < «K j D j (t), j - 1, 2 


(28) 


Let Dj - |»! - Xj)* 1 (|«j| _1 + D 2 - |t x - T 2 f l * Kf 1 )- 

■ 

Note that < a^. Thus if < 0 we readily see that is a steady 

state bound for D^(t). From these inequalities we derive, as In [2], the 
following conclusions: 

1. If Oj < 0 and jv^j _< M^, then the steady state solution of 

(11) satisfies the following inequality: 


,-l 


M. + 6K, D 


1“1 


h-U 

whereas If a 2 < 0 end |u p | _< M 2 , the steady state solution u m 
of (12) satisfies the Inequality 

|”-| < «2 + ® 2 D 2 

2. If » 0, the perturbation tern P^(t) can grow at most linearly. 

3. If dj > 0, the order of growth of P^(t) cannot exceed exp(a^t); 
note that the order or magnitude of all nonzero solutions of (19) or 
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(21) cannot exceed exp(a^t). 

Since oar assumptions Imply that f^(t) and f^(t) are bounded, we have 
therefore shown that the study of the boundedness of the solutions of (11) or 
(12) reduces to the study of the boundedness of the solutions of (19) or (21). 

If Oj < 0 and o 2 < 0 we shall say that the system (11), (12) (or 

(l)-(6)) is stable, if o^ - 0 and <» 2 < 0, or < 0 and <* 2 » 0, that the 

system has reached the stability boundary, and If o^ > 0 or > 0, that 
the system is unstable. Thus the system is stable if all its solutions are 
bounded. 


2.3. Estimates for and B 2 

tfe obtain estimates for the 6^ in terms of the coefficients of (11) and 
(12) and the signs of the a^. These estimates yield a simple method for 
deter min ing the stability of the system. Since (19) is identical with [2, 
(7)], we know the following: 

1/2 

1. If a x < 0, then Bj/Cj < Bj < (Aj + Kj) , and aj < 0. 


2. If o x - 0, then B /(^ 


3 , 


(A. + K.) 1/2 , and a! < 0. 


1/2 

3. If Oj > 0, then (Aj + 1^) 7 < Bj < B^Cj. 
If C 2 is real (i.e.. If o>I^ * 0), we also have: 


If 

°2 < °» 

then 

1/2 

®2 /C 2 < P 2 < (A 2 + K 2 ) » and 

a 2 < 

If 

° 2 • 0, 

then 

1/2 

®2 /C 2 “ P 2 " (A 2 +K 2 > * and 

“2< 

If 

« 2 > o. 

then 

(a 2 + k 2 ) 1/2 < b 2 < b 2 /c 2 . 



From these conclusions we also infer that if f ^ (t) and f 2 (t) are bounded, 
and Ij » 0, then the system (1) - (6) is stable if and only if 
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1/2 1/2 
B 1 /C 1 < (A x + Kj) , and B 2 /C 2 < (A 2 + 1^) . 

2.4 Resonance 

Proceeding as In [2], we readily see that If f ^ (t) - A^exp(lwt), 

j » 1, 2, then (1) - (6) can be In resonance only if a * 0 or 

1/2 

Q s /C s - [2a(bK s + aK b )/(I 1 + I 2 >] . 

3. Harmonic Analysis of the solutions. 

3.1. Preliminaries. 

We now study the properties of the Fourier transforms of the solutions. 

Following standard practice, we consider a tine interval of the form (c, d) 

0 < c < d < ». Let g^ C,d \t) » g(t) if c < t < d, and let g^ C,d \t) - 0 
otherwise. Thus, if F denotes the Fourier transform operator, we have: 
F[g^ C,d *](s) - (2v) ^ 2 £ g(t) exp(lst)dt, 
and proceeding as in [2] we see that 

lin F[v h (c,d) ](s) - 0, and lin Flu h (c,d) ](s) - 0, (29) 

and therefore 

lin F[v.^ C,d ^](s) • 0, and lin F[v 2 ^ C,d ^l(s) » 0, (30) 

e»-«» 1 c-** 

We want to study the properties of the graphs of the absolute values of 
p[ v (c» d )j(g) and F[u^ C * d ^l(s). From (29) it is clear that in order to 
obtain useful Information we need to study the Fourier transforms of the 
perturbation terns P^(t). 

3.2. Analysis of the perturbation terns. 

Let Qj(t) - f* exp[ Xj( t-x)]pj(x)dx, and 
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EL^t) - /* exp[ X^( fc-x) Jp j (x)dx, and let Q 2 (t), R 2 (t) be similarly defined in 
terms of Yj» y 2 » and p 2 (x). Clearly 

Pl (c,d) (t) - K 1 (x 1 -x 2 ) ” 1 [ Q J c,d)( t ) - aj c,d) (t)]. 

He have: 

F[Qj C,d ^l - (2x) 1>f2 exp[X 1 (t-x)lp 1 (x)dx exp(-sti)dx 

- (2x) 1 ^ 2 f d f* exp(-X.x)p. (x)dx exp[t(X. - si)]dt, 

' C ' O XI X 


and integrating by parts we obtain: 
F[Q< C ’ d) ](s) - 


(2x) 1 ^ 2 (X 1 -8i) 1 ^ 2 [exp(-dsi)p 1 (d) - expC-csDp^Cc)] - exp(-s(i)p 1 (t)dtj 
=» M 1 (S, c, d)/(A - si), 

where, since jp^t^ .1 

jMjts, c, d)| < (2 x )" 1/2 (2 + d - c)6. 

Using the same argument we also see that 
F[Rj C,d) ](s) - M^s, c, d)/X 2 - si). 


(31) 


where M^s, c, d) satisfies an inequality similar to (31). Thus, 
p[p (c,d)j( s ) _ K 1 (X 1 - X 2 )“ 1 [M 1 (s, c, d)/(X 1 -si) - M 2 (s, c, d)/(X 2 - si)]. (32) 


Similarly, 


F[P. 


.- 1 . 


2 C,d) ](s) - K^y^) X [M 3 ( 8 , c, d)/(y 1 -si) - M 4 (s, c, d)/(Y 2 -si)], 

v -1/2 


(33) 


where ^(s, c, d)| and |m 4 (s, c, d)| are bounded by (2ir) (2+d-c)6. 

Note, moreover, that there is no reason why M^(s, c, d) or M 3 (s, c » **) 
should vanish as c ♦ provided we keep the difference d-c constant. 

Since Vj = (l/2)(v-u) and v 2 ” (l/2)(v+u), for j ■ 1, 2 we obtain: 


V WV 


(c,d). 


where L Is a particular solution of (7) or (8), lim F[T 9 l(s) =* 0, 
j c ^ 
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(34) 


F[Hj C » d) ](s) - <l/2)(F[p< C ’ d) l(s) - F[p< C » d) ](a)), 

and 

F[H^ C,d) Ks) - (l/2)(F[p^ C,d) ](s) + F[P< C » d) ](s)). (35) 

Since X - si * + (0^ - si), and " si 3 + (B 2 - si), this aeans 

that the only value for which F[H^ C,d ^](s) and F[H^ c,d ^](s) may diverge as 

<Xj -*• 0 is s ■ Sj, and the only value for which they may diverge as 

«*2 0 is s 3 B 2 - 

Let o 3 max (a^, a^}. We identify the nonlinear natural frequency with 
those frequencies at which the PSD plots have relative maxima that become 
unbounded a -*• 0 . The method we have used to reach the conclusions of this 
section Is much simpler than the one we employed in [2]. 

3.3.3 Conclusions 

If o * max {a^, a^} and £ 3 B if 0 3 a or £ 3 B 2 otherwise, we 

conclude that as 0 ♦ 0 F[Vj C,d ^](s) and F[v 2 C,d ^](s) may diverge only at 

s 3 £ (or, if both -*■ 0 and o 2 ■*■ 0 simultaneously, at both B and 

B'), that for 0 negative and constant, but sufficiently close to zero, the 
gi-aphs of the absolute values of F[v^(c,d) J(s), j 3 1, 2 will have spikes 

near s 3 5 (or, in PSD plots, near C/sit), and that the magnitudes of these 

spikes need not decrease with time (l.e. , as c > «). 

4. Examples 

We now study the behavior of the solutions of (11), (12) for various 

2 

rotating speeds. We assume that 1^ 3 0, I 2 * 1.65 lbf.in.sec , 5 3 0.0015 in., 
2 

m 3 0.0587 lb.sec /in., a 3 4.15 in., b 3 7.1 in. If <f> denotes the angular 


XXXI I I- 10 



speed of the shaft In Hz., then K 3 0.0398<J> lb/ In, Q - 0.0201$ lb/ in, 

S 8 

C * 0.00001361b/ in, and K =» -1,800$ + 105,480,000. We also assume that 
s b 

_6 2 

e yl » 5.7 x 10 lb. sec /in, e y2 = -e yl , e z i “ e z 2 “ °» v (°) “ u (°) “ 6 » and 
v'(0) - u* (0) - 1. These values have been obtained from an actual bearing 
tester. Since, as we have already shown, when = 0 we have 
1/2 

Bj/Cj =* (Aj + K^) , we readily see that the value of $ for which this 

happens is <|> c ^ - 939 Hz. « 56,359 rpa, and that < 0 if $ < $ cl » Since 

1^ ■> 0, we can apply a similar procedure to conclude that 
$ c2 ■ 916 Hz. * 54,950 rpa, and that a 2 < 0 if $ < If f^ denotes 

the frequency that corresponds to the value of when = 0, it is 

readily seen that ■ 245 Hz. If f 2 denotes the frequency that 
corresponds to when a 2 ■ 0, we see that also <j> 2 = 245 Hz. Since a 2 

vanishes before a^, we reach the stability boundary when the shaft's rotating 
speed is 54,950 rpa. 

In Figs. 3 through 8 we see PSD plots for v 2 for various values of <)> 

ranging from 30,000 rpa to 57,000 rpa. (The plots for are siailar). To 

obtain these plots we first solved (11), (12) using a fourth order Runge-Kutta 
aethod. We then applied a Fast Fourier algor itha. The plots were obtained 
using 256 points and linear Interpolation, and are for the tine Interval 
5.120 sec. < t < 5.632 sec. The frequencies are measured in Hz. 

Since the aass Imbalance is so small, the forcing frequency $ is 
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undetectable. He see that as <fr Increases, the location of the nonlinear 
natural frequency a reaalns at 246 Hz. The magnitude of the spike Increases 
steadily, until around 50,000 rpa It starts to cl lab steeply. These exaaples 
show that the nonlinear natural frequency aay appear well before the stability 
boundary is reached. They also show that the location of this frequency Is 
not a good Indicator of stability aarglns, and that the approach to the 
stability boundary is accoapanled by a steep Increase In the size of the 
spike. These conclusions are similar to those we reached for a simple 
Jeffcott model In [2]. 
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